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Abstract 

Recent MEG results of a search for the lepton flavor violating (LFV) muon decay, /i — )■ 67, 
show 3 events as the best value for the number of signals in the maximally likelihood fit. 
Although this result is still far from the evidence/discovery in statistical point of view, 
it might be a sign of a certain new physics beyond the Standard Model. As has been 
well-known, supersymmetric (SUSY) models can generate the /i — )■ 67 decay rate within 
the search reach of the MEG experiment. A certain class of SUSY grand unified theory 
(GUT) models such as the minimal SUSY SO(IO) model (we call this class of models 
"predictive SO(IO) models") can unambiguously determine fermion Yukawa coupling ma- 
trices, in particular, the neutrino Dirac Yukawa matrix. Based on the universal boundary 
conditions for soft SUSY breaking parameters at the GUT scale, we calculate the rate 
of the /i — )■ 67 process by using the completely determined Dirac Yukawa matrix in two 
examples of predictive SO(IO) models. If we interpret the 3 events in MEG experiment 
as a positive signal and combine it with other experimental constraints such as the relic 
density of the neutralino dark matter and recent results on muon g — 2, we can pin down 
a parameter set of the universal boundary conditions. Then, we propose benchmark spar- 
ticle mass spectra for each predictive SO(IO) model, which will be tested at the Large 
Hadronic Collider. 
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The search for new physics beyond the Standard Model (SM) has been performed in 
a variety of energy scales. In addition to the direct searches for new physics at the Large 
Hadronic Collider (LHC), being the collider experiment with the highest energy at present, 
a search for the LFV processes at low energies is also very important in discovering new 
physics. This is because the LFV processes are highly suppressed in the SM, and any 
positive signal once observed can be an evidence of new physics. 

Recently MEG collaboration [T] has reported new results of a search for the /i — )■ 67 
decay and a maximally likelihood analysis sets an upper limit at 90% CL. on the branching 
ratio, BR(yU — )■ 67) < 1.5 x 10~^^, which is at the same level of the current smallest limit 
set by the MEGA experiment [2], BR(/i — )■ 67) < 1.2 x 10~^^. Very interestingly, the 
results of the MEG experiment also show 3 events as the best value for the number of 
signals in the maximally likelihood fit, which corresponds to pQ 

BR(/i ^ 67) = 3 X 10"^^ (1) 

for the center value. Although this result is still far from the evidence/discovery in 
statistical point of view, it might be a sign of a certain new physics beyond the Standard 
Model. 

Non-zero masses and fiavor mixings of neutrinos observed through the neutrino os- 
cillation phenomena imply that the lepton fiavor in each generation is not individually 
conserved, and therefore, LFV processes among the charged-leptons exist. However, in 
simply extended models of the SM so as to incorporate massive neutrinos, the rate of the 
LFV processes are strongly suppressed and far out of the reach of the experimental de- 
tection. This is because of the GIM mechanism which leads to a suppression factor of the 
ratio between the neutrino mass scale and the electroweak scale. It has been known [31 ll| 
that SUSY models can generate the experimentally detectable rate of the LFV processes 
through the LFV sources in soft SUSY breaking terms. 

The minimal supersymmetric Standard Model (MSSM) is one of the most promising 
candidate for new physics, providing a solution to the gauge hierarchy problem of the SM. 
In addition, the successful gauge coupling unification at Mqut — 2 x 10^^ GeV supported 
by the low energy data of the SM gauge couplings strongly suggests the paradigm of SUSY 
GUT. Thus, the SUSY GUT models with unified gauge groups are well-motivated theories 
at high energies. Among them, models based on the gauge group SO(IO) are probably the 
most compelling ones in terms of neutrino physics, because quarks and leptons in each 
generation are unified into a single 16 representation along with a right-handed neutrino, 
and the smallness of the neutrino masses can be naturally explained through the seesaw 
mechanism [5]. Furthermore, a class of SO(IO) models such as the minimal SUSY SO(IO) 
model [6] has a very interesting feature. Because of the complete unification of quarks 
and leptons into a single 16 representation and introduction of the minimal set of Higgs 
multiplets, fermion Yukawa matrices are highly constrained, and it is very non-trivial to 
fit all the current data of fermion masses and mixing angles including the neutrino sector. 
It has been shown [7] that the minimal SO(IO) model can simultaneously reproduce all the 
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observed quark-lepton mass matrix data involving the neutrino oscillation data. It is very 
interesting that after the data fitting, no free parameters are left and hence, all the fcrmion 
Yukawa matrices, in particular, the neutrino Dirac Yukawa matrix, are unambiguously 
determined. 

In this paper, we consider such a class of SUSY SO(IO) models (predictive SO(IO) 
models). We assume a suitable Higgs sector which breaks SO(IO) to the MSSM gauge 
group at the GUT scale with generating masses for the right-handed neutrinos. Below the 
GUT scale, the low energy effective theory of the model is described as the MSSM with 
the right-handed neutrino chiral multiplets. In this effective theory, the superpotential in 
the leptonic sector is given by 

Wy = + Y:^{e''^)djHa + \mr^.{v%Uv%)^ + (2) 

where the indices i, j run over three generations, i/„ and denote the up- type and down- 
type MSSM Higgs doublets, respectively, and M^.^. is the heavy right-handed Majorana 
neutrino mass matrix. We work in the basis where the charged-lepton Yukawa matrix 
Ye and the mass matrix Mr^. are real-positive and diagonal matrices: Y^^ = YaSij and 
Mji-. = diag{Mji^, Mr2, Mr^). Thus the off-diagonal components of the neutrino Dirac 
Yukawa coupling matrix Y^, break the conservation of the lepton flavor in each generation. 
Note that once we fix a predictive SO(IO) model, all components in the neutrino Dirac 
Yukawa matrix Y^, and the mass spectrum of the right-handed neutrinos arc completely 
determined. The soft SUSY breaking terms in the leptonic sector is described as 



- ^soft = ^1 (^1) ij + ^Ri i^l) ^Rj + ^Ri (^g) ,j ^Rj 



+ {A'Ai^jH^ + ^iei^jHa + h.c) 

+ (]-MiBB + ]-M2W''W'' + ^MaG^G" + h.c^ . (3) 

Current experimental results severely constrain soft SUSY breaking parameters to be 
almost flavor blind and real. In our analysis, we adopt the most commonly considered 
scenario, the constraint MSSM (CMSSM), and impose the universal boundary conditions 
for soft SUSY breaking parameters at the GUT scale: 

Ml = Ma = Ma = M1/2 , (4) 

and extrapolate the soft SUSY breaking parameters to the electroweak scale according 
to their renormalization group equations (RGEs). Although the boundary conditions 
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are flavor-blind, the LFV terms in the soft SUSY breaking parameters such as the off- 
diagonal components of (jnj^ and A^J are induced by the LFV interactions, namely, the 

neutrino Dirac Yukawa couplings. For example, the LFV effect most directly emerges in 
the left-handed slepton mass matrix through the RGEs such as 




where the first term in the right hand side denotes the normal MSSM term with no LFV. 
In the leading- logarithmic approximation, the off-diagonal components {i ^ j) of the 
left-handed slepton mass matrix are estimated as 

where the distinct thresholds of the right-handed Majorana neutrinos are taken into ac- 
count by the matrix L = \og[MGUT / Mji-]6ij . 

The effective Lagrangian relevant for the LFV processes {ii — )■ £-,-7) is described as 

^eff = -l^iM.^F^" {A^Pl + A^Pr) , (7) 

where Prl = (1 ± 75)/2 is the chirality projection operator, and A^ r are the photon- 
penguin couplings of 1-loop diagrams in which chargino-sneutrino and neutralino-charged 
slepton are running. The explicit formulas of Al^r etc. used in our analysis are summa- 
rized in [8]. The rate of the LFV decay of charged-leptons is given by 

m ^ ^.7) = ^nii (|Aif + ) . (8) 

The following approximation formula [8] is useful to understand the order of magnitude 
of the LFV decay rate: 

-=2 /ri' \ 2 I i^^f) P 

r(..^.,,)._™f,x(^) ^-^W^, (9) 

where Ms is the average slepton mass at the electroweak scale, and (^Am|^ is the slepton 

mass estimated in Eq. ([6]). From these formulas, once we obtain the information of the 
neutrino Dirac Yukawa coupling matrix and the right-handed neutrino masses, we can 
predict the LFV decay rate, for a given set of universal boundary conditions and tan /3. 

Recent various cosmological observations, in particular, the Wilkinson Microwave 
Anisotropy Probe (WMAP) satellite [9], have established the ACDM cosmological model 
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with a great accuracy. The rehc abundance of the cold dark matter (CDM) in 2a range 
has been measured as 

f]^^j^^/i2 = 0.1131 ±0.0034. (10) 

As is well-known, in SUSY models with the i?-parity conservation, a neutralino, if it 
is the lightest sparticle (LSP), is a promising candidate for the CDM in the present 
universe. When we apply the WMAP result to the relic density of the neutralino LSP, 
the parameter space of the CMSSM is dramatically reduced into the narrow stripe due 
to the great accuracy of the WMAP data [TU]. We take into account this cosmological 
constraint in our analysis of the LFV processes below. For tan /3 = 45, /i > and Aq = 0, 
for example, we can find the approximate relation between mo and M1/2 such as 

along which the observed relic abundance, QcoAih"^ = 0.113, is realized. In our analysis, 
we have employed the SOFTSUSY 3.1.4 package [T^ to solve the MSSM RGBs and 
produce mass spectrum. Then, the relic abundance of the neutralino LSP is calculated 
by using the micrOMEGAs 2.4 [12] with the output of SOFTSUSY in the SLHA format 

Now let us consider two example models in the class of predictive SO (10) models. The 
first example is the minimal SUSY SO(IO) model analyzed in [7], in which only one 10 
and one 126 Higgs multiplets have Yukawa couplings with 16 matter multiplets such as 

Wy = Yi^l6iH,ol6j + Yli,16,HuGlQj, (12) 

where 16^ is the matter multiplet of the i-th generation, Hiq and H12Q are the Higgs 
multiplet of 10 and 126 representations under SO(IO), respectively. By virtue of the 
gauge symmetry, the Yukawa couplings, Yio and Y12Q, are, in general, complex symmetric 
3x3 matrices. 

We assume some appropriate Higgs multiplets, whose vacuum expectation values 
(VEVs) correctly break the SO(IO) gauge symmetry into the Standard Model one at 
the GUT scale, Mqut — 2 x 10^^ GeV. Note that Hiq and H12Q, respectively, includes one 
pair of up-type and down-type Higgs doublets in the same representation as the pair in 
the MSSM. Suppose that the Higgs doublets in the MSSM are realized as one light linear 
combination of the two pairs of Higgs doublets and the other combination is as heavy 
as the GUT scale, realizing the MSSM as the low energy effective theory. The Higgs 
doublet in the representation (10,1,3) under the Pati-Salam subgroup, G'422 =SU(4)cX 
SU(2)iX SU(2)/j, is embedded in the 126 representation, whose nonzero VEV breaks the 
Pati-Salam subgroup to the SM one and at the same time, generates the masses of the 
right-handed neutrinos. Providing suitable VEVs for the Higgs multiplets, the quark and 
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lepton mass matrices are characterized by only two basic mass matrices, Mio and M126, 
and three complex coefficients, Cio, Ci26 and cr, at the GUT scale: 

= CioMio + Ci26Mi26 

Md = M10 + M126 

Md = CloMio - 3Ci26Mi26 
Me = Mio-3Mi26 

Mr = CrMuQ, (13) 

where M^, M^, M^, Me, and Mr denote the up-type quark, down-type quark, neutrino 
Dirac, charged-lepton, and right-handed neutrino Majorana mass matrices, respectively. 
Except for cr, which is used to determine the overall neutrino mass scale, this system has 
fourteen free parameters in total [Hj, and the strong predictability to the fermion mass 
matrices. 

In Ref. [7], thirteen electroweak data of six quark masses, three mixing angles and 
one phase in the Cabibbo-Kobayashi-Maskawa (CKM) matrix, and three charged-lepton 
masses are extrapolated to the GUT scale according to the RGEs with a given tan/3, and 
are used as inputs at the GUT scale. Solving the GUT mass matrix relation among quarks 
and charged- leptons obtained by Eq. f|T3|l . we can describe the neutrino Dirac mass matrix 
Md and M126 as functions of only one free parameter a, the phase of the combination 
(3io + Ci26)/(— Cio + Cue). It has been shown [7J that the appropriate value of a and 
Cr can reproduce the light neutrino mass matrix consistent with the observed neutrino 
oscillation data (as of 2002). For tan/3 = 45 and cr = 3.198 fixed, the right-handed 
Majorana neutrino mass eigenvalues are found to be (in GeV) [7] 

Mr^ = 1.64 X 10^\ Mr^ = 2.50 X 10^^ Mr^ = 8.22 x 10^^ (14) 

where cr is fixed so that Am^ = 2 x 10~^eV^. In the basis where both of the charged- 
lepton and right-handed Majorana neutrino mass matrices are diagonal with real and 
positive eigenvalues, the neutrino Dirac Yukawa coupling matrix at the GUT scale is 
unambiguously determined and explicitly given by 

/ -0.000135 - 0.00273i 0.00113 + 0.0136i 0.0339 + 0.0580i \ 
y;= 0.00759 + 0.0119i -0.0270 - 0.00419i -0.272 - 0.175i . (15) 
V -0.0280 + 0.00397i 0.0635 - 0.0119z 0.491 - 0.526z / 

Although we will use the data in Eqs. ([Hj) and f|T5|) for our numerical analysis below, 
the neutrino oscillation parameters given by these equations are, unfortunately, more 
than 3a away from the current neutrino oscillation data [12]. Some extension of the 
minimal SO(IO) model is necessary to improve the data ffiting for the neutrino oscillation 
parameters [TB]. However, the resultant neutrino Dirac Yukawa couplings and the right- 
handed neutrino mass spectrum are not so much changed by this improvement, and we 
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believe that our result for the LFV process in this paper gives, at least, the correct order 
of magnitude for the prediction of a slightly extended minimal SUSY SO(IO) model. 

As another example of predictive SO(IO) models, we consider a simple SUSY SO(IO) 
GUT in five dimensions (5D) proposed in [17], which can ameliorate several theoretical 
problems in the minimal (renormalizable) SUSY SO(IO) model (see [17] for detailed dis- 
cussions on the theoretical problems of the minimal SUSY SO(IO) model). This model 
is defined in five dimensions with the fifth dimension compactified on the S^jiZ^ x Zg) 
orbifold [18]. The SO(IO) gauge symmetry in five dimensions is broken by the orbifold 
boundary conditions to the Pati-Salam (PS) symmetry SU(4)cX SU(2)ixSU(2)/j. All 
matter and Higgs multiplets are arranged to reside only on a brane (PS brane) at one 
orbifold fixed point where the PS gauge symmetry is manifest, so that low energy effective 
four dimensional description of this model is nothing but the PS model with a special set 
of matter and Higgs multiplets. The PS symmetry is broken at the normal GUT scale, 
Mgut — 2 X 10^^ GeV, in usual four dimensional manner, while the full gauge coupling 
unification is realized after incorporating threshold corrections of Kaluza-Klein modes 
in the bulk gauge multiplets. Phenomenology of sparticles [19] and applications of the 
model to cosmology such as inflation [20] and baryogenesis via leptogenesis [21] have been 
investigated. 

In this model with a simple set of Higgs multiplets, the structure of fermion Yukawa 
couplings are almost the same as the one in the minimal SUSY SO(IO) model. Assuming 
appropriate VEVs for Higgs multiplets, fermion mass matrices are obtained as (see [17] 
for the detailed definition of the mass matrices and the coefficients) 

= CioMi,2,2 + Ci5Mi5,2,2 , 

Md = Mi,2,2 + Mi5,2,2 , 

Md = CioMi,2,2 - 3Ci5Mi5,2,2 , 

Me = Mi,2,2 - 3Mi5.2,2 , 

Mr = Ci?Mio,l,3, (16) 

which are characterized by three fundamental mass matrices Mi^2,2, Afi5^2,2 and Mio,i,3 and 
complex coefficients cio, C15 and cr. Comparing these expressions with Eq. f lT3|) . we can 
see that the combination of two mass matrices of Mi 2,2 and M15 2,2 among M^, M^, Mjy, 
and Me is the same as that of Miq and M126 in the minimal SUSY SO (10) model and, 
therefore, the procedure for fitting the realistic Dirac fermion mass matrices is the same 
as in the minimal SO(IO) model. On the other hand, Mr is fully independent on the 
above four Dirac fermion mass matrices, whereas in the minimal SUSY SO(IO) model 
it is described by M126 and not independent. This fact enables us to improve the data 
fitting of the neutrino oscillation parameters. 

Through the seesaw mechanism [S], the light neutrino mass matrix is given by 

= YjMj^Vu = UmnsDMLns (17) 
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in the basis where the mass matrix of charged lepton is diagonal with real and positive 
eigenvalues. Here f„ is the VEV of the up-type Higgs doublet in the MSSM, D^, is the 
diagonal mass matrix of light neutrinos, and Umns is neutrino mixing matrix. This is 
equivalent to the expression of the right-handed neutrino mass matrix as 



Mn = vl {YMTBAjD^'U^BMyu 



Once the information of the Dirac Yukawa coupling, the mass spectrum of the light neu- 
trinos, and the neutrino mixing matrix is obtained, we can fix the right-handed neutrino 
mass matrix. In order to determine M^, we follow the manner in [2T] . 

We consider the normal hierarchical case for the light neutrino mass spectrum, for 
simplicity, and describe Di, in terms of the lightest mass eigenvalue mi and the mass 
squared differences: 

= diag (mi, y^Amfa + mj, \J Amj^, + mf ^ . (19) 

Here we adopted the neutrino oscillation data [T5] : 

Amf^ = 7.59 X 10"^ eV^ Amj^ = 2.43 x 10"=^ eV^ (20) 

In addition, we assume the neutrino mixing matrix of the so-called tri-bimaximal form 
[22] 



/ /I \ 

V 3 V 3 ^ 



TBM 



1 /I /I 
3 V 2 
1 A 



\ Ve V3 \J 2 ' 



(21) 



which is in very good agreement with the current best fit values of the neutrino oscillation 
data |15j. As has been discussed above, the data fit for the realistic charged fermion mass 
matrices is the same as in the minimal SO(IO) model, and we here use the numerical 
value of Yy in Eq. (IT^ at the GUT scale for tan /3 = 45, in the basis where the charged 
lepton mass matrix is diagonal. Substituting Eqs. ( TT^ . ( TT^ . ( 12U]) and (1^ into Eq. (ITSil . 
we obtain the right-handed neutrino mass matrix as a function of only mi. Since it has 
been shown [21] that the simple 5D SO(IO) model can reproduce the observed baryon 
asymmetry of the present universe for mi ~ 1.8 x 10^'^ eV through non-thermal leptoge- 
nesis, we take this value as a reference value for mi. In this way, is now completely 
determined, but not yet diagonalized. Changing the basis to diagonalize M/j, we find the 
neutrino Dirac Yukawa matrix, 

-0.00119 + 0.0000268i -0.00108 - 0.000485i -0.000392 + 0.000421i 
Yy=\ 0.00135 + 0.00167i -0.0253 + 0.00154i 0.0237 + 0.00085^ | , (22) 
-0.0265 - 0.0173i 0.0609 + 0.0275i 0.790 - 0.0436i 



7 



and the diagonal right-handed neutrino mass matrix with eigenvalues (in GeV) 

Mr^ = 1.03 X 10^°, Mr^ = 7.55 X 10^\ Mr^ = 3.22 x 10^^ (23) 

Now we are ready to analyze the /i — )■ 67 decay rate by using the completely determined 
neutrino Dirac Yukawa coupling matrix and the right-handed neutrino mass eigenvalue^, 
Eqs. ([T3D and f[T^ for the minimal SO(IO) model while Eqs. fl22]) and fl2^ for the simple 
5D SO(IO) model. For tan /3 = 45, Aq = and /i > 0, the branching ratio of yU — )■ 67 for 
the minimal SO(IO) model is shown in Fig. 1 as a function of the universal gaugino mass 
at the GUT scale, along the relation of Eq. (ITTl) to satisfy the observed relic density for 
the neutralino dark matter. The short dashed line corresponds to the branching ratio of 
Eq. ([T]) obtained from the MEG results if the three events are considered as a positive 
signal, while the long-dashed line to the upper limit of the branching ratio set by the 
MEGA experiments. The universal gaugino mass M1/2 — 790 GeV (which corresponds 
to niQ ~ 415 GeV) reproduces the MEG results. 

Fig. 2 depicts the same result as in Fig. 1, but for the simple 5D SO(IO) model. In this 
case, the universal gaugino mass reproducing BR(/i — > 67) = 3 x 10~^^ is relatively light, 
Ml/2 — 420 GeV (which corresponds to mo — 272 GeV). This is because the components 
in the neutrino Dirac Yukawa matrix relevant to the LEV between 1st and 2nd generations 
are smaller than those in the minimal SO(IO) model and a lighter sparticle mass spectrum 
is necessary in order to achieve the same branching ratio (see Eq. ([9j)). 

For each model, we have pinned down the soft SUSY breaking mass parameters in the 
CMSSM, in the light of the MEG results. The set of the CMSSM parameters proposes a 
good benchmark point for the SUSY search at the LHC, and we present (sparticle) mass 
spectra for each model in Table 1, along with other observables which can be compared 
with the current experimental bounds. In both results, the lower bound on the Higgs 
boson mass rrih > ll^A GeV [21] is satisfied. Other phenomenological constraints we 
consider here are 

2.85 X 10-^ < BR(6 ^ 57) < 4.24 X 10"^ (2(t) [25], (24) 

BR{Bs fi^fx-) < 5.8 X 10-8 [26]. (25) 

The results for the minimal SO(IO) model satisfy these constraints, while BR(6 — )■ 57) for 
the simple 5D SO(IO) model is marginal (about SAa away from the center value). 

The muon anomalous magnetic dipole moment (muon aMDM) has been measured in 
a great precision as [27] 

aJ^^P = 11659208.0(6.3) x 10"^°, (26) 

where the number in parentheses shows la uncertainty. On the other hand, the SM 
predictions were calculated [28j (see also [29 | [30l [3T| [32l [33] ). 

a^^[T] = 11659193.2(5.2) x 10"^°, 
^ See |23j for previous analysis on the minimal SO(IO) model with more general parameter sets. 
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a? [e^e-] = 11659177.7(5.1) x 10"^°, (27) 

by using data in the hadronic r decay and e+e~ annihilation to hadron, respectively, in 
calculating the hadronic contributions to the muon aMDM, aSM[,-] 

and a^^[e"'"e ]. The 
deviations of the SM predictions from the experimental result are given by 

Aa^Jr] = P - a^^lr] = 14.8(8.2) x 10"^°, 
Aa^[e+e-] = af^'P - af^[e+e-] = 30.3(8.1) x 10"^°, (28) 

which correspond to 1.8cr and 3.7cr deviations, respectively. These deviations may be from 
the sparticle contributions. Table 1 also include the sparticle contributions to Aa^. We 
can see that the result for the minimal 5D SO(IO) model favors the deviation obtained by 
using the data of the hadronic r decay, while the deviation from the e~^e~ data is favored 
by the result for the simple 5D SO(IO) model. 

With Ml/2 (and mo) determined in each SO(IO) model, we can also predict the LFV 
decay rates of tau lepton in the same analysis as for fj, — >■ 67. For the minimal SO (10) 
model, we find 

BR(r ^ /i7) ~ 5.74 x 10"^° and BR(r ^ 67) ~ 1.32 x 10"^°, (29) 

while for the simple 5D SO (10) model 

BR(r ^ /i7) ~ 5.53 x 10"^° and BR(r ej) ~ 1.34 x 10"^°. (30) 

These predicted values are well below the current experimental upper bounds [M], BR(r — i- 
A^7) ^ 4.4 X 10~^ and BR(r — )■ ej) < 3.3 x 10~^. Note that improvements by one or two 
orders of magnitude are expected at future super B-factories p5j and the branching ratio 
of order 10"^'' can be tested in the near future. 

In summary, recent MEG results of a search for the /i — )■ 67 decay might show a sign 
of a certain new physics beyond the Standard Model. The primary candidate for new 
physics which can give the rate of the LFV decay accessible to the MEG experiment is 
supersymmetric models. A certain class of SUSY SO(IO) GUT models (predictive SO(IO) 
models) can unambiguously determine fermion Yukawa coupling matrices, in particular, 
the neutrino Dirac Yukawa matrix and the mass spectrum of the right-handed neutrinos. 
As an example of such models, we have considered the minimal SO(IO) model and a 
simple 5D SO(IO) model. The effective theory of the SO(IO) models below the GUT scale 
is described as the MSSM with the right-handed neutrinos. In these models, even if the 
universal boundary conditions for the soft SUSY breaking parameters at the GUT scale, 
LFV terms are induced in the sparticle sector through the RGE effects with the neutrino 
Dirac Yukawa matrix, which give rise to a sizable rate of the — )■ 67 decay. Since the 
neutrino Dirac Yukawa matrix and the mass spectrum of the right-handed neutrinos are 
completely determined, the LFV rate can be predicted once the set of input parameters 
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in the CMSSM is given. Requiring to reproduce BR(yU — > 67) = 3 x 10~^^ given by the 
MEG results and combining this requirement with the cosmological constraint on the 
observed rehc abundance of neutrahno dark matter, we have pinned down the parameter 
set in the CMSSM and then, have proposed benchmark points for each SO(IO) modeL 
We have checked other phenomeno logical constraints for each benchmark point. Finally, 
our strategy in this paper is applicable also to other SUSY (GUT) models if they can 
completely determine the neutrino Dirac Yukawa matrix and the right-handed neutrino 
mass spectrum. Depending on models, different benchmark points can be pinned down, 
which will be tested at the LHC. 
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Figure 1: The branching ratio BR(yU — )■ 67) for the minimal SO(IO) model as a function 
of Ml/2 (GeV) along the cosmological constraint of Eq. f ill I) . The short dashed line 
corresponds to the MEG result, BR(/i — )■ ej) = 3 x 10~^^ when three evens are considered 
as a positive signal, while the long dashed hne to the upper limit on BR(/i — )■ ej) = 
1.2 X 10"^^ given by the MEGA experiment. 
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Figure 2: The same figure as Fig. 1 but for the simple 5D SO(IO) model. 



13 



Table 1: Mass spectra (in GeV) and phenomenological constraints for the two SUSY 
SO (10) models with the universal boundary conditions in the CMSSM 





minimal SO (10) model 


simple 5D SO (10) model 


mo 


415 


272 


Ml/2 


790 


420 


^0 








tan/3 


45 


45 


ho 


119 


115 


Ho 


786 


449 


Ao 


787 


449 




791 


457 


9 


1756 


981 


Al,2,3,4 


333,631,928,938 


171,324,535,548 


Xl,2 


631, 938 


324, 548 




1576, 1645 


898, 934 




1582, 1643 


901,931 


fcl,2 


1409, 1473 


784, 849 


ii,2 


1266, 1475 


698, 864 




667,667,619 


386, 386, 355 




511,672 


317, 395 




342, 642 


186,392 


BR{b S7) 


3.27 X 10"^ 


2.36 X 10"^ 




1.04 X 10~^ 


4.95 X 10-9 


Aa^ 


12.0 X 10-1° 


37.7 X 10-1° 




0.113 
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